Lensing Time Delays and Cosmological Complementarity 
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Time delays in strong gravitational lensing systems possess significant complementarity with dis- 
tance measurements to determine the dark energy equation of state, as well as the matter density 
and Hubble constant. Time delays are most useful when observations permit detailed lens modeling 
and variability studies, requiring high resolution imaging, long time monitoring, and rapid cadence. 
We quantify the constraints possible between a sample of 150 such time delay lenses and a near term 
supernova program, such as might become available from an Antarctic telescope such as KDUST 
and the Dark Energy Survey. Adding time delay data to supernovae plus cosmic microwave back- 
ground information can improve the dark energy figure of merit by almost a factor 5 and determine 
the matter density fim to 0.004, Hubble constant h to 0.7%, and dark energy equation of state time 
variation Wa to 0.26, systematics permitting. 



I. INTRODUCTION 

Complementarity between cosmological probes in- 
creases their leverage on the cosmological model parame- 
ters, crosschecks results through differing systematic un- 
certainties, and breaks degeneracies. These all play im- 
portant roles in elucidating the nature of our universe: 
the energy densities in matter and dark energy, the scale 
of the universe through the Hubble constant, and the 
characteristics of the dark energy behind the current cos- 
mic acceleration. 

Most probes, however, have substantial similarity in 
their parameter dependencies, involving the same combi- 
nations of ingredients entering into the Hubble parameter 
as a function of redshift, H{z). Distances (and volumes) 
in particular, are essentially equivalent, and growth of 
structure also depends similarly on H{z). Looking for a 
high degree of complementarity, especially to determine 
the dark energy equation of state value and time varia- 
tion, [ij investigated the use of distance ratios present in 
strong gravitational lensing as a means of breaking this 
degeneracy. 

While lensing distance ratios involve the mass struc- 
ture of the lens, and so are not purely geometric, there 
has been impressive progress in modeling the mass distri- 
butions in lensing systems (e.g. and so it is worth 
considering strong lensing distances in more detail as a 
cosmological probe, in particular for its complementarity. 
Here we revisit [l| with several important distinctions: 1) 
we concentrate on time delays, due to the recent obser- 
vational successes d, [H and modeling advances; 2) we 
consider a more realistic range of future observational 
prospects, involving projects starting to get underway, 
which will have important implications for complemen- 
tarity; and 3) we carry out studies of the science reach 
as a function of redshift range, and in the presence of 
spatial curvature. 

Several authors (e.g. [6l-[l2j) have addressed the sta- 
tistical power of strong lensing time delays from further 
future surveys such as LSST, calculating the numbers of 
lenses found and with measured time delays, and project- 



ing possible Hubble constant or cosmological constraints. 
These, however, treat the lensing systems as an ensemble 
to average over, and in fact identify the mass modeling 
as a major uncertainty capable of degrading constraints 
substantially. Here we concentrate on what are some- 
times called "golden lenses" , although now the meaning is 
not systems with some special symmetry but rather ones 
where the survey design has specifically provided data 
enabling detailed construction of the lens mass model. 
The number of such systems will be much less but we 
find they can have significant scientific leverage. 

In Section [IT] we discuss the cosmological impact of 
time delay measurements and their complementarity 
with other probes. Section Hill considers reasonable pos- 
sibilities for survey data sets in terms of number and red- 
shift range of time delay systems, and analyzes their con- 
straints in conjunction with a mid term supernova survey 
and cosmic microwave background data. Survey require- 
ment issues with respect to imaging resolution, time sam- 
pling, etc. for time delay measurement, lens modeling, 
and systematic error control are outlined in Section IIVI 
with specific reference to the Antarctic optical/infrared 
telescope program. 



II. TIME DELAYS AS COSMIC PROBE 

Strong gravitational lensing causes multiple images of 
distant sources, with the light rays from the images tak- 
ing different amounts of time to propagate to the ob- 
server. The time delay involves two parts: a geometric 
delay from different path lengths and a gravitational time 
delay from traversing different values of the gravitational 
potential of the lens. Thus both the image positions and 
the lens mass model must be accounted for. Time de- 
lays are observed by looking for coordinated variations 
in the flux from the images, e.g. of time varying quasars, 
which requires long time, well sampled monitoring. Typ- 
ical galaxy lens induced delays are ~ 60 days and the 
desired measurement accuracy is a couple of days or bet- 
ter. 
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To translate the image angular positions into spatial 
positions, for computing both the path length and the 
gravitational potential effects, one needs the (conformal) 
distance to the lens, ri, to the source, r^, and between 
the source and lens, ris. Only in flat space is r/^ = rg—ri. 
The particular combination of distances central to time 
delays is 

T=^. (f) 

ris 

Specifically, following Q, the time delay of an image 
at position 9 on the sky relative to an unlensed source at 
position j3 is 

M{eJ) = ''-^{i + zi)cj,{eJ), (2) 

ris 

where the distance ratio is T, containing the key cosmo- 
logical information, and </> is the Fermat or time delay 
potential given by 

^0J)=^±lf^-^^O). (3) 

We see that the first term on the right hand side is the 
geometric delay and the second term is the lensing poten- 
tial delay with V^tA ~ 2k for k the dimensionless lensing 
projected surface mass density. 

The time delay potential (p connecting T to the ob- 
served time delays therefore depends on the lens mass 
distribution, the model for which is built up from infor- 
mation on image positions and flux ratios and perhaps 
surface brightness morphologies, ideally from many im- 
ages including arcs. See d, foi" details on the mod- 
eling process and calculation of the potential factor. All 
the cosmological information comes from T (cf. the ap- 
proach of [l3|), with the uncertainties in the potential 
factor entering into the error propagation, together with 
measurement uncertainties. 

The time delay probe T has interesting properties with 
regards to the cosmological parameter leverages. As 
noted by [l|, the sensitivities to dark energy parameters 
wq and Wa , where the dark energy equation of state is well 
fit by w{z) — wq + Waz/{1 + z) [14|, are actually posi- 
tively correlated in contrast to standard distance mea- 
surements such as from Type la supernovae. This offers 
hope for complementarity with such probes. Further- 
more, the sensitivity to the dimensionless matter den- 
sity rim at low redshift is remarkably low compared to 
solo distances, leading to the possibility of breaking the 
usual degeneracy between matter density and dark en- 
ergy equation of state. Finally, the ratio depends linearly 
on the Hubble scale Hq ^ , and since [l^ researchers have 
sought to use lensing time delays to measure the Hubble 
constant. 

In Figure [T] we highlight these special properties. The 
derivatives dhiT/dp give the sensitivities for each pa- 
rameter p and are exactly what enters into a Fisher ma- 
trix analysis for cosmological parameter estimation. Raw 



(unmarginalized) sensitivities can be read directly, e.g. 
dhiT/dp/O.Ql = 10 means that a 1% measurement of T 
delivers an uncertainty a{p) = 0.1. This must be folded 
in with the covariances between parameters: sensitivity 
curves having the same (reflected) shape indicate highly 
anticorrelated (correlated) parameters. For the time de- 
lay probe, the curve shapes are not very similar - a good 
sign for breaking degeneracies - and we see the unusual 
positive correlation between wq and Wa over the whole 
range zi = 0-0.6 (where for simplicity we have assumed 
Zs = 2zi). 
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FIG. 1. Sensitivity of the time delay distance combination 
T = rirs/ris to the cosmological parameters is plotted vs 
lens redshift. Curves with opposite signs at the same redshift 
indicate positive correlations between those parameters - very 
unusual for the dark energy equation of state variables wo and 

Wa- 



To take advantage of the odd correlation properties 
of T to give strong complementarity in probing cosmol- 
ogy, we include Type la supernova distances as another 
probe, having very different degeneracies. The super- 
nova distance-redshift relation has no sensitivity to the 
Hubble constant h = iJo/(100km/s/Mpc) however, this 
being convolved with the unknown supernova absolute 
luminosity. To proflt from the time delay dependence on 
/i, therefore, we also use cosmic microwave background 
(CMB) information, which determines the physical mat- 
ter density combination ilmTi^ very well (but not par- 
ticularly Vlrn by itself). This further offsets the weak 
dependence of T alone on the matter density, and so the 
weakness of each is turned into strength in complemen- 
tarity. 
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III. COSMOLOGICAL LEVERAGE 

Another interesting property of the time delay probe 
is that its useful and unusual eorrelation properties oc- 
cur at low redshift, for zi = 0-0.6. Detailed observa- 
tions of lensing systems will be easier there, where the 
lens galaxy and source images will not be as faint as at 
higher redshift. We therefore take as our baseline a sur- 
vey producing time delay measurements at zi = 0.1-0.6 
(there is relatively little volume for lensing systems below 
zi = 0.1), and then study variations of this. For simplic- 
ity we fix Zs ~ 2zi; although there will be a distribution 
of source redshifts this has little impact on the cosmology 
estimation (see, e.g., §5.5 of 0) and we have exphcitly 
checked that using instead Zg = Azi affects the dark en- 
ergy figure of merit (uncertainty area) result by less than 
1%. In most of this section we assume a spatially flat 
universe, studying the effect of an additional parameter 
for curvature in Section [III Bl 



A. Cosmological Parameter Constraints 

To the time delay measurements we add supernova dis- 
tance (SN) and CMB information and carry out a Fisher 
matrix analysis to estimate the cosmological parameter 
constraints. For the supernovae, we take a mid term 
sample reasonable for the next five years, consisting of 
150 SN at z = 0.03-0.1 from the Nearby Supernova Fac- 
tory [l^, 100 SN per 0.1 bin in redshift from z = 0.1-1 
as from the Dark Energy Survey (DES: [l3|) with fol- 
low up spectroscopy, and 42 SN between z = 1-1.7 as 
from Hubble Space Telescope observations such as the 
CLASH [3 and CANDELS ^ surveys. This seems like 
a reasonable estimate for a mid term, well characterized 
supernova sample. Each supernova is given a 0.15 mag 
(7% in distance) statistical uncertainty and each redshift 
bin of 0.1 has a systematic floor at dwsys = 0.02 (1 -I- z) 
added in quadrature to the statistical error. Thus the 
supernova sample is systematics limited out to z = 1. 
For CMB data, we take Planck quality information con- 
sisting of determination of the geometric shift parame- 
ter R to 0.2% and the physical matter density ^mh'^ to 
0.9%, roughly corresponding to constraints from the lo- 
cation and amplitude, respectively, of the temperature 
power spectrum acoustic peaks. The parameter set is 
{Q.m,wo,Wa,h^ AA.}, where ^A is the convolution of the 
supernova absolute luminosity and the Hubble constant. 

Current measurements can deliver the time delay 
probe r to 5% for a lensing system, dominated by 
systematic uncertainties for individual systems. With a 
survey designed to find many strong lensing images and 
characterize them accurately, it may be possible to con- 
sider 1% measurements of T in each redshift bin of 0.1 
from z = 0.1-0.6. This can be thought of as either 25 
strong lenses per bin (150 total), or fewer lenses with bet- 
ter accuracy than 5% per system from a survey designed 
to gather data needed to control systematics, or a com- 



bination of the two. We discuss the survey requirements 
in Section HVl 

Figure [2] shows the dramatic improvement in the dark 
energy equation of state parameters (marginalized over 
the other parameters) when adding the time delay probes 
of 1% accuracy over z = 0.1-0.6. The area of the er- 
ror contour in WQ-Wa tightens by a factor 4.8 over that 
from SN+CMB alone. All the cosmological parameters 
are better determined by factors of 2.6-3.1. Time delays 
therefore have great complementarity with the supernova 
and CMB probes, and such a strong lensing survey would 
be highly valuable scientifically. 
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FIG. 2. 68% confidence level constraints on the dark en- 
ergy equation of state parameters wq and Wa using mid term 
supernova distances and CMB information, and with (solid 
curve) or without (dashed curve) time delay measurements. 
The time delay probe demonstrates strong complementarity, 
tightening the area of uncertainty by a factor 4.8. 

The absolute level of the constraints with time delays is 
impressive as well. The Hubble constant is determined to 
0.0051, or 0.7%; the matter density to 0.0044 (1.6%), 
and the present value of the dark energy equation of state 
Wo to 0.077 and its time variation Wa to 0.26. While 
falling short of the results from a space survey of su- 
pernovae (with CMB), such a mid term program could 
deliver important insights into the nature of cosmic ac- 
celeration and the cosmological model. 

The baseline time delay sample adopted seems plausi- 
ble, but let us consider variations to see how the cosmo- 
logical constraints depend on the survey characteristics. 
It may be difficult to find enough strong lens systems 
at the lowest redshifts, due to the limited volume. Note 
however that the SLACS survey has been successful in 
detecting lenses (20j . if not necessarily measuring time 
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delays, at zi « 0.1, and this depends on the source pop- 
ulation targeted. Nevertheless, if we cut the time de- 
lay information to the range z = 0.3-0.6 (so 100 strong 
lenses), we find that this reduces the figure of merit (in- 
verse uncertainty area) by 25%. The greatest effect is 0.5 
on the Hubble constant determination, since this is what 
low redshift time delays excel at, with a{h) degrading by 
55%. This then propagates into the constraint, which 
weakens by 41%. These can be somewhat ameliorated if 
we have some information from z = 0.1-0.2, e.g. using 12 
rather than 50 time delays in this range recovers almost 
half the constraining power. 

Conversely, suppose that a strong lensing time delay 
survey could be extended out to z/ = 1, still at the 1% 
accuracy per 0.1 redshift bin. Then the figure of merit 
improves by 40%, though the constraints on fi^ and h -0.5 
only gain by 6%. Detailed characterization of the lens- 
ing systems at such high redshift could be problematic, 
however, due to lower fluxes and signal to noise. The red- 
shifts of well characterized time delay systems is slowly 
being pushed out toward z/ = 1 [2l|, [2^ . 
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B. Including Curvature 

Spatial curvature enters together with the Hubble pa- 
rameter into cither the angular or luminosity distance 
between observer and source. Degeneracy between the 
curvature density H.^ = 1 — ilm — ^de and dark energy 
equation of state can be severe; for example see Fig. 6 of 
[231 for effects on wq, Wa or for general w{z). This can 
be broken by using a wide redshift range of distances; in 
particular high and low redshift distance measurements 
can separate the curvature density from other compo- 
nents. Another possibility is direct measurement of the 
Hubble parameter as well as distances (e.g. from the ra- 
dial baryon acoustic oscillation scale), or distance ratios 
appearing in gravitational lenses or large scale structure 
(see, e.g., |25l - j27| ). This has the advantage of not neces- 
sarily requiring high redshift measurements. 

We now examine the role that time delay measure- 
ments can play in breaking the curvature degeneracy, if 
the universe is not assumed to be spatially fiat. Figure [3] 
shows the results when we allow for curvature in the cos- 
mology fitting, using time delays, supernovae distances, 
and CMB information. 

The dark energy equation of state uncertainty indeed 
degrades, with the area figure of merit declining by a fac- 
tor 4.1. This shows the degeneracy is not fully broken, 
but should be contrasted with the factor 20.2 degrada- 
tion from including curvature with only the SN+CMB 
data for constraint. Thus the time delay probe is a use- 
ful tool even/especially when allowing for spatial curva- 
ture. Most of the covariance affects the time variation 
Wa, with its uncertainty doubling. The present value wq 
is only determined 12% worse, and the errors on fl^ and 
h increase by 29% and 27%. The curvature itself is esti- 
mated to (T(Ofe) = 0.0063. 



FIG. 3. 68% confidence level constraints on the dark energy 
equation of state parameters wo and Wa using time delay, 
mid term supernovae, and CMB information, assuming spa- 
tial flatness (solid curve) or allowing curvature (dashed curve). 
The time delay probe moderates the curvature degeneracy, 
restricting the area degradation to a factor 4, rather tlian 20 
without time delay data. 



IV. SURVEY CHARACTERISTICS 

In order to use time delays as a cosmological probe in 
the individual lensing system approach, the survey must 
deliver detections and accurate measurements of the time 
delays, detailed modeling of the lens systems, and con- 
trol of other systematic uncertainties. Systematics in- 
clude microlensing that induces variability, differentially 
altering the images' light curves, and projected mass not 
truly part of the lens, altering the mass modeling. 

To detect a large sample of time delay systems, a wide 
field survey is needed, but to characterize them through 
accurate image positions, splittings, and flux variations 
requires high resolution imaging. Interestingly, a tele- 
scope at an excellent seeing site such as Dome A, Antarc- 
tica [2^ could fulfill both roles. The Kunlun Dark Uni- 
verse Telescope (KDUST: [2§|), a 2.5 meter telescope 
planned for Dome A would be situated above the low 
ground layer and possibly have 0.3" median seeing in the 
optical, 0.2" in the low background noise infrared. 

The advantages of high resolution imaging for strong 
lensing are crucial and manifold [9, 30]. Such seeing also 
helps to separate the images from contamination by the 
lens galaxy light. To take advantage of this excellent see- 
ing for strong lensing, the point spread function (PSF) 
would need to be stable, or algorithms developed to fit 
simultaneously the PSF and lens mass model. The stable 
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winter weather, with low winds and large isoplanatic an- 
gle, at Dome A could be advantageous. KDUST surveys 
would overlap with Dark Energy Survey fields, as well as 
those of the South Pole Telescope and LSST. DES could 
supply much of the supernova sample, although super- 
nova programs, at either low or very high redshift, are 
also being studied for KDUST [3l|. 

Measuring time delays accurately from detected strong 
lensing systems requires a long time baseline, since the 
time delay distribution of interest is in the range of ~10~ 
100 days. The long Antarctic night offers advantages here 
of continuity, although end effects from the long Antarc- 
tic day mean that not all systems at the upper end of this 
distribution will be usable, in particular the longer clus- 
ter lenses (which also likely have larger external mass ef- 
fects). Dense time sampling enables accurate time delay 
determination, and ameliorates the effect of microlensing 
systematics, and again the Antarctic site allows continu- 
ous viewing of fields and regular sampling, every 24 hours 
or even more often. (Indeed, detection of time delay per- 
turbations can be used as a probe of dark matter sub- 
structure and properties [32j.) While photometric red- 
shifts are likely good enough to remove most projected 
mass contamination, follow up spectroscopy for accurate 
determination of the redshifts of the images and lens con- 
stituents is necessary. A spectrographic telescope is being 
considered for Dome A, but spectroscopy is needed in any 
case for DES fields (e.g. for the supernovae), so KDUST 
gains further advantage from synergy with DES. 

The prospects for 1% measurement of the time delay- 
redshift relation in a 0.1 bin in redshift for lenses at 
z = 0.1-0.6 seem reasonable. Improvements in control of 
systematics could tighten the current 5% accuracy, and 
such surveys will build the statistics as well. The num- 
ber of time delay systems baselined in this article - 150, 
a one order of magnitude increase over current levels - is 
plausible, as is the range of supernovae data, making the 
science case for time delay surveys of interest for further 
investigation. 

V. CONCLUSIONS 

We have quantified the significant complementarity 
as cosmological probes that strong gravitational lensing 
time delays, involving distance ratios, have with solo dis- 
tance measurements such as from Type la supernovae. 
A well designed time delay survey can add to practical, 
near term supernova and CMB data to provide surpris- 
ingly incisive constraints on the dark energy equation of 
state, the Hubble constant, and the matter density. The 
improvement in equation of state area uncertainty (figure 
of merit) is almost a factor 5 over the data sets without 
time delays. 

Time delays also significantly ameliorate the degenera- 
cies in parameter determination caused by allowing for 
spatial curvature, again improving the area uncertainty 
by a factor 5. Determination of the Hubble constant to 



0.7% as well would be valuable for several astrophysical 
and cosmological applications. 

We have focused on what seem to be near term, rea- 
sonable data sets. An exciting possibility for achieving 
these is telescopes being developed at promising Antarc- 
tic astronomical sites, such as KDUST at Dome A. If 
these truly deliver high resolution, stable seeing much 
better than conventional ground based optical conditions 
(if not quite space quality) , the baseline time delay survey 
considered here to deliver one order of magnitude times 
larger sample of well characterized time delay systems ap- 
pears practical. Another advantage is the synergy with 
other southern surveys, such as the Dark Energy Sur- 
vey in the near term. (While we have intentionally not 
extrapolated to long term developments, synergy with 
LSST is clear as well.) 

Systematic uncertainties would be ameliorated by the 
high resolution imaging, whether single epoch to char- 
acterize in detail the lens model and separate the host 
galaxy light, or multiepoch to finely measure the flux 
variations and measure clean and accurate time delays. 
The redshift range for the survey could be modest, 
zi K, 0.1-0.6, and we presented how the cosmological con- 
straints change if a narrower or wider range is considered. 

In attempting to stay within straightforward practical- 
ity we have not discussed exciting ideas such as testing 
for deviations from general relativity. After all, the same 
principle is used in the solar system, with spacecraft sig- 
nal time delays providing stringent limits on other grav- 
ity theories; also some screening mechanisms that re- 
store gravity to general relativity are expected to kick 
in on scales accessible to cosmological strong lensing (ssl - 
[35| . Many other astrophysical applications exist for a 
high resolution imaging survey (especially with low noise 
in the infrared), such as using strong lenses as gravita- 
tional telescopes to study early structures and the epoch 
of rcionization. 

Complementarity between cosmological probes offers 
the strongest and most robust leverage for revealing the 
scale and contents of our universe, and the nature of the 
cosmic acceleration. The combination of time delays, su- 
pernova distances, and CMB data provides an exciting 
level of insight with near term surveys. 
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